Disorder induced roughening transition of many elastic lines in a periodic potential 
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The competing effect of a periodic pinning potential and random point disorder is studied for arrays 
of elastic lines or directed polymers. The groundstates are investigated by exact combinatorial 
optimization. In both two and three dimensions a phase diagram is found with two or three distinct 
phases: a strictly flat phase if the disorder is bounded and weak, a weakly fluctuating phase for 
intermediate valleys depths, where the lines roughen individually on a scale smaller than the line-line 
distance, and a rough phase for strong disorder, where the roughness follows the scaling with pure 
point disorder. The line wandering in the transverse direction leads in the three-dimensional, rough 
phase to an entangled state with a complicated topology. 



Ensembles of elastic lines form one of the most inter- 
esting examples of the interplay of disorder and an or- 
dering tendency. Much of the attention is motivated by 
the experimental connections to dirty type-II supercon- 
ductors, in which flux lines (FL) interact with point or 
columnar defects ||l|-3|. This interaction gives often rise 
to beneficial, technologically important effects since the 
impurities pin lines, individually or collectively and cre- 
ating them in a controlled fashion has been demonstrated 
by various techniques. The theory of arrays of FLs with 
a wide variety of possible disorder backgrounds has re- 
ceived lots of interest Q. It is worth noting that the 
physics of an individual FL has connections to paradig- 
matic questions in non-equilibrium statistical mechanics 
[H,|j through mappings to the Kardar-Parisi-Zhang equa- 
tion of kinetic roughening in surfaces, and to the Burgers' 
equation of vortex- free turbulence H . 

Here we analyze with exact numerical tools the rough- 
ening or disordering of arrays of elastic lines at zero 
temperature when point disorder competes with a pe- 
riodic potential that tries to order the array into a reg- 
ular structure. The FL analogy is an Abrikosov lattice 
in the presence of point disorder. For both two- and 
three-dimensional (2D and 3D) arrays of lines we find 
a transitions in terms of the ratio q = (e) /A, where A 
measures the depth of the potential valleys, and (e) is 
the average strength of the point disorder. The picture 
we gain from our study is in contrast to recent renor- 
malization group studies of a model for elastic periodic 
media with point disorder in the continuum limit pO| , pl| . 
Based on these one would expect for an arbitrary dis- 
order strength q in two dimensions asymptotically (in 
system size) a rough state with diverging FL displace- 
ments and line-to-line displacement correlations, whereas 
we find this to be true only above qc2 ■ In other words, the 
lattice model for a flux line array with a periodic poten- 
tial shows a roughening transition which should, at least 
in 2d, be absent according to analytical theories for an 
apparently related continuum model. This is also surpris- 
ing since for the same system without periodic potential 



good agreement between analytical [[l2| and numerical 
p3| results has been achieved. After presenting the nu- 
merics we will discuss the apparent contradiction with 
analytical results. 

We are interested in an assembly of N elastic lines 
described by the Hamiltonian 

^ - E / dz{im' + E,(#.) ^nt[r.(z) - r,(z)] 
+ K[r,(z),z]+yp[r,(z)]| . 



(1) 



Yi{z) S TV^^"^ (with d — 2 and d — 3 studied here) is 
a displament vector, while z is the longitudinal coordi- 
nate in a system of height H; K[r, z] describes the point 
disorder, which we take to be delta-correlated with vari- 
ance e; Vint [r — r'] is a short-range repulsive interaction 
between the lines (e.g. hard-core) and V^[r] a periodic 
potential with period a in all transverse space directions. 
We assume its minima to be well localized, i.e. they have 
a width that is small against the interaction range of the 
lines. This allows only single occupancy of the potential 
valleys and we concentrate on the case in which the line 
density (p — L'^^^ /N) is such that each potential valley 
is occupied by exactly one elastic line, i.e. p = l/a'*~^. 
A similar model has been studied in d = 1 + 1 in 111] 
for the competition between point disorder and random 
columnar defects, i.e. a potential Vp that is not periodic 
but random. 

We study a lattice version of the continuum model (|l|) 
that lives on the bonds of a square lattice [15|. The struc- 
ture of the potential is depicted in Figurej^. We have 
valleys of effective depth —A, and bulk disorder from a 
probability distribution P{eij) for the bonds of the lattice 
model. P is taken to be either a flat, bounded distribu- 
tion so that (e) — e,„aj,/2, or an unbounded exponential 
distribution with (e) = e, the decay constant of the dis- 
tribution. Each bond can only be occupied by a segment 
of a single line (hard core interactions) and the occu- 
pation of each bond (ij) costs a particular amount of 
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FIG. 1. Periodic potential in 2D. The depth of the val- 
leys is denoted by A and the nearest neighbour distance by 
a. Additional point disorder e accomplishes the energy land- 
scape. The FL can only enter and leave the system via the 
energetically neutral arcs connecting the source respectively 
the sink with the potential valleys. 
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j + Ay where the Ay 's are zero inside the 
potential wells and constant elsewhere. This reproduces 
also the elastic energy, since all bonds cost some positive 
energy. The lattice Hamiltonian then reads 



iJ(x) 



(y) 



(2) 



where X^fji") i^ ^ sum over all bonds (ij) joining site i and 
j of a d-dimensional, e.g. rectangular {L'^~^ x H) lattice, 
with open boundary conditions (b.c.) in all space direc- 
tions. The FL enter the system via the plane z = Q and 
leave it via the plane z = H, see Fig. |l|. The FL configu- 
ration is defined by the set of variables Xij that can take 
on the values Xij — 1 (if a line is passing bond (ij) and 
Xij = 0). For the configuration to form lines on each site 
of the lattice the configuration x has to be divergence 
free (i.e. V • x = 0, where V- denotes the lattice diver- 
gence). All top-sites are attached to an extra site, via 
energetically neutral arcs, e = 0: the source with diver- 
gence V -x = +N, where N is the number of FLs, and all 
bottom sites to the target t with V • x = —N. Finding 
the minimum energy configuration of N lines now be- 
comes equivalent to a minimum cost flow problem with 
the energy function (0) and the mass balance constraints 
expressed in the site divergences. This task can be solved 
exactly in polynomial time (complexity 0{N ■L'^~^H) by 
applying the successive-shortest-path algorithm which is 
described in detail in ||l5|-[l7| . 

Figure g demonstrates with a series of snapshots the 
geometry involved in the calculations, and the typical 
behavior with increasing g in 2d and 3d. In both cases 
the lines are pinned to the energetically favorable valleys 
for small q, and finally for large q a cross-over to a rough 
state takes place. In 3d one can observe that the lines 
wander almost freely in such conditions. 

We discriminate between the different regions in the 




FIG. 2. Optimal ground state configurations in 2D (top) 
and 3D (bottom) for different point disorder strengths q, in- 
creasing from left to right. In the flat phase (left) the FL are 
trapped completely inside the potential valleys. 



phase diagram by looking at the behavior of the average 
transverse fluctuation or roughness w of the lines: 



w{L,H) 






, (3) 



where Tj ~ H^^ L dzYi{z) and [. . .]av denotes the dis- 
dorder average. By studying very large longitudinal sys- 
tem sizes H > 10^ we are able to extract the saturation 
roughness w{L) = lim^f^oo w{L, H) for a finite system 
of transverse size L. Note that we have chosen open 
b.c: the transverse fluctuations cannot get larger than 
the system size. Other quatities of interest are the size In 
of the longitudinal excursions (the average distance be- 
tween the locations at which a line leaves a valley and 
returns to it); and the total number of potential valleys 
PV that a line visits between its entry and terminal point 
in the limit if — > c». 

In Fig. |3|we show our data for the roughness w and ^|| 
as a function of 1/g in 2d, in Fig. |j w and PV as a func- 
tion of q in 3d; both for bounded disorder. In Fig. |5| we 
show our data for w for unbounded disorder. The picture 
that emerges is the following. In the flat region we have 
w{L) = 0, Z|| = and PV — 1, i.e. the lines lie completely 
in the potential valleys. This region q < qd exist only 
for bounded disorder. For the uniform distribution no en- 
ergetically favourable transverse fluctuation can exist as 
long as g < A. That qd > 1 follows from the fact that 
we are at full occupancy, N — Ny where Ny is the num- 
ber of valleys, for <z < <Zci ~ 2 the groundstate consists 
always of N straight lines regardless of dimension. For 
unbounded disorder this flat region does not exist, since 
the probability for a sequence of high-energy bonds in 
the valleys that pushes the lines out of it is always pos- 
itive. In the weakly fluctuating region for gci ^ <? ^ qc2 
the Hues roughen locally. Here one has w > and Z|| > 0, 
independent of the systems size L, and PV = 1. The 
transverse fluctuations of flux lines are bounded by the 




FIG. 3. Roughness w in 2D as a function of disorder 
strength q for bounded disorder, gci and qc2 are shown. In 
the flat phase w = 0, whereas w > for q > qd . No transver- 
sal system size dependence is observered in the weakly fluktu- 
ating phase. The inset shows Z||. Each data point is averaged 
over n = 20 (L = 128) up to n = 600 (i = 8) disorder 
configurations, a = 4. 



average line distance or valley separation a. The cen- 
tral feature is that lines fluctuate individually, so that a 
columnar defect competes with point disorder. Both in 
2D and in 3D a strong columnar pin strictly localizes the 
line [|l reducing the line-to-line interaction to zero. 

With increasing q the transverse fluctuations increase, 
and when their lengthscale l± becomes comparable to 
the inter-line distance a the physics changes. When in- 
dividual lines can jump from one valley to a neighboring 
valley, i.e. exactly when PV > 1, a collective rearrange- 
ment of the whole line ensemble becomes possible and 
manifests itself in a system size dependence of the satu- 
ration roughness (see Fig. Isl for 2d and H for 3d) . Thus 
the rough phase is characterized by a dependence of w 
and ^j[ on the transverse system size and by PV > 1. 
In particular the latter criterion facilitates the numerical 
determination of the location of the roughness transition 
qc2- Wc flnd that the critical disorder strength qc2{a) 
increases monotonically with a. 

In the limit g ^ oo the collective behavior of the 
transverse Auctions of the lines crosses over to that of 
a line array in the absence of the periodic potential, 
though the scaling of displacement correlations and the 
system roughness is slightly decreased compared to the 
case without the potential because of the lengthscale of 
the pinned parts of the lines (the remaining tendency to 
localize in the valleys). For an elastic medium model 
this would imply that the correlations should depend in 
a logarithmic fashion on distance. In 2d, in particular, it 
should be w '^ In L |lj,|l^,|l^ and our data indicate that 
this scaling is indeed obtained above qc2 ■ 




FIG. 4. (Top) The same as in Fig. |, now in 3d. (Bottom) 
The FL start to jump from one potential valley to the next 
at qc2 so that the number of potential valleys each line visits 
changes from pu = 1 to p« > 1 at this threshold (bottom). 
n = 6 {L = 128) up to n = 100 (L = 5), a = 4. 



A fundamental difference exists between the real-space 
geometry of two-dimensional and three-dimensional sys- 
tems. In 3D and for small line densities, the lines can 
wander freely since there is no hard-core expulsion in 
contrast to 2D. This implies that in the thermodynamic 
limit, the line ensemble becomes entangled. We char- 
acterize this by computing the fraction of lines N^nt/N 
that, when forced to start and end in an arbitrary one of 
all the potential minima, enter and exit in different val- 
leys. In the limits H -^ cxd, L —> oo Nent{q)/N develops 
a flrst-order jump at the roughening transition from zero 
to unity. N^nt = N implies that all the lines are entan- 
gled, though in a finite system the saturation value of 
Nent/N is limited by the finite probability that at least 
one out of N integers ends up at exactly at the same po- 
sition in a random permutation. This finite size feature 
can be seen in the data presented in Fig. |6[ 

Finally we relate our findings to periodic elastic media 
with point disorder ||ll| . These are supposed to represent 
a coarse grained description of flux line arrays such as d3) 
and predict the absence of a roughening transition in 2d 
contrary to our results. The reason for this discrepancy 
is that such models operate in a different limit of many 
lines between the potential valleys and interactions that 
renormalize to an elastic constant for the line displace- 
ments ll^. Simple scaling considerations show [O that 
this model in d=l+l is always rough, essentially due to 
the fact that the energy for terrace excitations of length 
scale L in the disorder free case (A — 0) has no bulk 
contributions. On the other hand, in the model (|l|) we 
have considered here excitations, e.g. steps of individual 
lines from one potential valley to the next, do have a 
bulk contribution, due to the full occupancy of the nar- 
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FIG. 5. Same as Fig. W for unbounded disorder. The 
inset shows 3-D. Though hardly visible from this data qd 
vanishes. But we still find a finite value of qc2, below which 
the roughness shows no L-dependence both, in 2D and 3-D. 
n = 10 up to n = 200, both in 2D and 3-D. 



row potential valleys. Therefore we expect a roughening 
transition for p) at a non- vanishing disorder strength, as 
confirmed by our numerical results in d—2,3 

To conclude, we have analyzed with the aid of exact 
combinatorial optimization methods ensembles of elas- 
tic lines (or directed polymers or FLs) in the presence 
of a confining periodic potential and competing random 
point disorder. The main finding is a transition between 
'rough' and 'flat regimes, in both two and three dimen- 
sions, at a finite potential strength. It arises since even 
rare fluctuations are not able to induce line-line inter- 
actions when the filling factor of the system is at unity. 
Introducing dislocations (in the form of missing lines) or 
extra lines is likely to complicate the physics, as in the 
rough substrate model. In the rough phase the physics is 
characterized by correlations that increase with system 
size. In three dimensions we find an entangled phase, in 
which the lines form a topologically complicated geomet- 
ric configuration. It needs to be characterized properly, 
and the effect on the glassy physics should be studied. 
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FIG. 6. Number of entangled lines N^nt normalized by the 
total number of lines. The horrizontal lines indicate the ana- 
lytically expected saturation values for random permutations 
of N lines. 
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